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Abstract 

Coronaviruses (CoVs) represent a global public health threat, exemplified by the severe acute respiratory 
syndrome (SARS) and Middle East respiratory syndrome (MERS) outbreaks. Using fecal samples collected 
from five bat species between 2014 and 2016 in Finland and RT-PCR, RT-qPCR, and NGS, we identified CoVs 
in 10 of 79 (13%) samples, including two novel bat species-CoV relationships. Phylogenetic analysis revealed 
Alphacoronavirus and Betacoronavirus species clustered among previously identified bat and human viruses. 
These results expand the known northern distribution and host species range of bat-borne CoVs. 
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Introduction 

T he emergence of severe acute respiratory syndrome 
(SARS) and Middle East respiratory syndrome (MERS) 
established coronaviruses (CoVs) as a major global public 
health threat (Drexler et al. 2014). Spillover from wildlife 
sources is implicated for these emergence events, with subse¬ 
quent research emphasis placed on documenting the wildlife 
host range and geographical distribution of CoVs to better un¬ 
derstand the origins of human disease, and help design spill¬ 
over mitigation strategies and improve public health readiness. 

CoVs are classified into four genera, designated Alpha-, 
Beta-, Gamma-, and Deltacoronavirus (Woo et al. 2012, 
Drexler et al. 2014). Gamma- and Deltacoronavirus are pre¬ 
dominately found in birds, whereas Alpha- and Betacor¬ 
onavirus infect mammals, including humans (Woo et al. 
2012). An especially high diversity of Alpha- and Betacor¬ 
onavirus has been identified in bat species, supporting the 
relatively recent recognition of bats as major wildlife reser¬ 
voirs for zoonotic viruses (Calisher et al. 2006) and indicating 
their ancestral host role in CoV evolution (Drexler et al. 2014). 
Despite this, minimal research has been conducted in northern 
ecosystems, including Fennoscandia (Finland and Scandina¬ 
via), where CoVs of bats have not been previously reported. 

Materials and Methods 

Noninvasive cross-sectional surveys were conducted in 17 
separate sites in southcentral Finland between 2014 and 


2016. Bats were captured using mist nets, harp traps, or by 
hand directly from roost sites, and individual fecal samples 
were collected from five bat species ( Eptesicus nilssonii, 
Myotis brandtii, Myotis daubentonii, Pipistrellus nathusii, 
and Plecotus auritus). Bats were kept in cotton bags for up to 
20 min and fecal pellets were collected from the bag or during 
examination. Bat species were identified using morphology, 
which is a sufficient method in Finland with only 13 known 
bat species. In addition to individual sampling, guano sam¬ 
ples were collected from known roosting sites. 

Samples were stored in RNAlater (QIAGEN, Hilden, 
Germany) in the field, and processed in two different labo¬ 
ratories using slightly different methods. RNA was extracted 
from the 2013 and 2014 samples using the QIAamp Viral 
RNA Mini Kit (QIAGEN) and primary screening was per¬ 
formed using the OneStep RT-PCR Kit (QIAGEN). We used 
previously described primers targeting the CoV RdRP-gene 
(11-FW 5' TGATGATGSNGTTGTNTGYTAYAA 3' and 
13-RV GCATWGTRTGYTGNGARCARAATTC 3'; Mur- 
adrasoli et al. 2009). Positive samples were confirmed using 
the same RT-PCR kit and different previously described RT- 
PCR primers with slight modifications ( Deltacoronavirus 
conserved primers, 5' GGTTGGGACTATCCTAAGTGTG 
3', and 5' CATCATCAGATAGAATCATCAT 3'; Woo et al. 
2006), which yields a 440 bp amplicon. 

After the project moved to a different laboratory with 
different reagents and established methods in 2015, RNA 
was extracted from the 2015 and 2016 samples using TriPure 
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Table 1. Total Number of Each Sample Type 
Collected from Each Bat Species in Finland (Number 
of Positive Samples/Total Number of Samples) 


Bat species 

Fecal 

samples 

Roost 

guano 

samples 

Total 

(prevalence 

%) 

Eptesicus nilssonii 

2/22 

0/14 

2/36 (6) 

Myotis brandtii 

2/6 

— 

2/6 (33) 

Myotis daubentonii 

4/27 

2/5 

6/32 (19) 

Pipistreilus nathusii 

0/1 

— 

0/1 (0) 

Plecotus auritus 

0/2 

0/1 

0/3 (0) 

Unidentified species 

— 

0/1 

0/1 (0) 

Total (%) 

8/58 (14%) 

2/21 (10%) 

10/79 (13%) 


(Roche. Basel, Switzerland) and screened using the qScript 
One-Step SYBR Green qRT-PCR Kit, Low ROX (Quanta 
Biosciences, Beverly, MA) based on the same primers 
used for the 2013 and 2014 samples (same 11-FW as 
mentioned and slightly modified 13-RV GCATWGTRT- 
GYTGNGARCARAATT 3', no probes; Muradrasoli et al. 
2009). Positive samples were reverse transcribed to cDNA 
using the RevertAid Premium Kit (ThermoFisher Scientific, 
Waltham, MA) and confirmed through a previously de¬ 
scribed nested RT-PCR method (de Souza Luna et al. 2007), 
which yields ~400 bp amplicon on the nested round, using 
a Phusion Flash High-Fidelity PCR Master mix (Thermo- 
Fisher Scientific). PCR amplicons were sequenced using 
Sanger sequencing. 



HQ184056.1 Myotis daubentonii, Spain 

EU375875.1 Myotis daubentonii, Germany 

MG923567.2 Myotis daubentonii, Lieto, 2014 

MG923574.1 Myotis daubentonii, Mustasaari, 2016 

MG923569.1 Myotis daubentonii, Kirkkonummi. 2014 

MG923573.2 Myotis daubentonii, Helsinki, 2016 

MG923568.2 Myotis daubentonii, Vihti, 2014 

MG923571.2 Myotis daubentonii, Kirkkonummi, 2014 

DQ648858.1 Scotophilus kuhlii. China 

EU375869.1 Pipistreilus nathusii. Germany 

MN065811.1 Myotis brandtii. Helsinki. 2016 

HM368166.1 Myotis myotis, Germany 

KY009625.1 Myotis ricketti. China 

KY009623.1 Myotis fimbriatus, China 

DQ249224.1 HKU6 Myotis ricketti, China 

AY567487.2 NL63 Human 

KY073744.1 Triaenops afer, Kenya 

KY073745.1 Triaenops afer. Kenya 

AF3044601 229E Human 

JQ410000.1 Alpaca repiratory coronavirus 

KY502384.1 Myotis emarginatus. Luxembourg 

MG923572.2 Myotis brandtii. Espoo, 2015 

EF544563.1 Myotis occultus. USA 

KF500940.1 Hypsugo savii, Italy 

KF500946.1 Pipistreilus kuhlii. Italy 

HQ184059.1 Hypsugo savii. Spain 

KY009626.1 Eptesicus serotinus, China 

KC869678.4 Neoromiaa ca pens is, South Africa 

JX869059.2 MERS Human 

KY009614.1 Myotis pequinius. China 

KY009634.1 Eptesicus serotinus, China 

HQ184062.1 Eptesicus isabellmus. Spain 

MG923570.1 Eptesicus nilssonii, Joroinen, 2014 

KF312399.1 Eptesicus serotinus. Italy 

DQ074652.1 HKU4 Tylonycteris pachypus, China 

DQ249217.1 HKU5 Pipistreilus abramus. China 

KC545386.1 European hedgehog. Germany 

NC_006577.2 HKU1 Human 

FJ647225.1 MHV Mouse 

AY391777.1 OC43 Human 

DQ011855.1 PHEV Pig 

NC_004718.3 SARS Human 

DQ022305.2 HKU3-1 Rhinolophus sinicus, China 

KY502395.1 Rhinolophus ferrumequinum. Luxembourg 

GQ404795.1 Rhinolophus hipposideros, Slovenia 


Alpha-CoV 


Beta-CoV 


FIG. 1 . Phylogenetic consensus tree indicating the evolutionary relationships among partial coding sequences of the CoV 
RdRP-gene (392 nucleotides) from Finnish bats and CoVs found elsewhere in bats and other species. Terminal node 
information includes the bat species, and the location and year of sample collection (year of only the samples reported in 
this study). Posterior probabilities of internal nodes are located next to given nodes and the scale bar shows expected 
nucleotide changes per site. Boldface nodes represent sequences reported in this study. The tree was made using MrBayes 
3.2.7ax 64 with a GTR + G nucleotide substitution model, two Markov chain Monte Carlo (MCMC) runs with four chains, 
1,000,000 samples, a sampling frequency of 1000, and a burn-in of 25%. CoV, coronavirus. 
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Extracted RNA from CoV-positive samples was also used 
for sequencing on a MiSeq platform (Illumina, San Diego, 
CA). Preparation of the samples and assemblies of the reads 
was done as described previously (Forbes et al. 2019). 
Quality of the assemblies was assessed by aligning reads to 
contigs using BWA alignment algorithm in UGENE (1.29.0 
64-bit version; Okonechnikov et al. 2012). The same method 
with reference sequences instead of contigs was used to re¬ 
trieve cytochrome b and cytochrome c oxidase subunit I se¬ 
quences from MiSeq reads of guano samples for host species 
identification (available upon request). 

Sequences acquired from Sanger and MiSeq sequencing 
were aligned (Muscle 3.8.31) with sequences of CoV RdRP 
from GenBank, after which the best substitution model was 
selected (MEGA7 7.0.26 Model Selection [ML] feature with 
default settings). A phylogenetic tree was constructed using a 
Bayesian framework with MrBayes (3.2.6x64; Ronquist 
et al. 2012). 

Results 

Ten of 79 total samples contained CoV RNA (13%) and at 
least 1 positive sample was obtained from three of the five 
different bats species, with CoV prevalence ranging from 
6% to 33% (Table 1). The exceptions were P. nathusii and 
P. auritus. 

CoV sequences were obtained from nine samples: two from 
M. daubentonii roost sites (both alpha-CoV, GenBank accession 
numbers MG923567.2 [776 bp] and MG923569.1 [395 bp]), 
four from M. daubentonii (all alpha-CoV, GenBank acces¬ 
sion numbers MG923568.2 [622bp], MG923571.2 [4805 bp], 
MG923573.2 [15 837 bp], and MG923574.2 [28 045 bp]), two 
new CoV-host relationships; M. brandtii (two alpha-CoVs, 
GenBank accession numbers MG923572.2 [436 bp], and 
MN065811.1 [28 119 bp]) and E. nilssonii (beta-CoV, Gen¬ 
Bank accession number MG923570.1 [395 bp]). The con¬ 
structed phylogenetic tree showed sequences belonging to 
Alphacoronavirus and Betacoronavirus genera, which were 
clustered among known bat and human viruses (Fig. 1). 

Discussion 

This report increases knowledge on the geographical dis¬ 
tribution and host species range of bat-borne CoVs. We show 
for the first time that bat-borne CoVs occur in northern 
Europe, and identified two previously unknown bat species- 
CoV associations. 

The alpha- CoV sequences from M. daubentonii , one of 
which was from a bat captured in Mustasaari (GenBank ac¬ 
cession number MG923574.2) on the west coast of Finland 
and another from a guano sample collected from bat colony in 
Lieto (MG923567.2) in the southeast corner of Finland, 
clustered among sequences from central and southern Eur¬ 
ope, supporting previous findings of gene flow among Eu¬ 
ropean populations of this host (Laine et al. 2013). Other CoV 
sequences from M. daubentonii were from the southern coast 
and formed their own cluster. 

Meanwhile, one of the alpha- CoVs from M. brandtii 
(MN065811.1) clustered among other CoVs of Myotis spp. 
from Europe. There is a possibility that this is due to a cross¬ 
species transmission, since another alpha- CoV from M. 
brandtii (MG923572.2) was phylogenetically separated from 
those sequences, with one exception being a sequence from 
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M. emarginatus from Luxembourg, which parallels phylo¬ 
genetic reconstructions made of the Myotis genus (Ruedi 
et al. 2013). The Beta- CoV from E. nilssonii clustered with 
Betacoronaviruses of other Eptesicus spp. Most of these 
findings indicate host fidelity and the codiversification of 
CoVs with their host species. 

Although we show that bat species in northern Europe are 
host to CoVs, potential human health concerns are prema¬ 
ture based on these initial findings, and must be balanced 
against the important environmental services that bat species 
provide. Bats play a critical role in ecosystem functioning, 
including insect pest management, seed dispersal, and pol¬ 
lination (Kunz et al. 2011), and are protected in Finland and 
most other parts of Europe. 

Our identification of bat-borne CoVs in northern Europe, 
including two novel host-CoVs relationships, will enable 
more in-depth study about the pathogens carried hy bats 
living in urban, rural, and forest habitats in northern latitudes, 
their evolutionary relationships with host species, and po¬ 
tential risks to human health. Further research, including 
whole genome targeting methods (e.g., NGS) from fresh 
samples from wider geographical and temporal ranges, is 
required to better understand the evolutionary history and 
epidemiological properties of these CoV-host relationships, 
and ultimately the risks they may pose to human health. 
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